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We have investigated skyrmions in iV-component quantum Hall systems. We find that SU(iV) 
skyrmions are the lowest energy charged excitations for filling factors v = 1, 2, . . . , N — 1 for small 
enough symmetry breaking terms. N > 2 skyrmions can be realized in Si QH systems based on 
the (110) or (111) interfaces of Si, or perhaps in Si (100) systems, where the spin and valley isospin 
together provide an SU(4)-symmetry, or in multilayer QH systems. We also present Hartree-Fock 
results for a phenomenological easy-axis SU(2)-breaking model appropriate to valley degeneracy. 
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Textured quasiparticle excitations in the quantum Hall 
(QH) effect involving nontrivial configurations of electron 
spin ("skyrmions") were proposed in jjj. There is now 
strong experimental evidence that such excitations are 
the lowest energy charged quasiparticles at filling frac- 
tion v = 1 (2|-Q. The physics behind these excitations 
is as follows. First, even in the absence of Zeeman cou- 
pling, the QHE system is an exchange ferromagnet, ow- 
ing to the degeneracy of the Landau levels. Second, spin 
textures carry charge; the charge density is proportional 
to the skyrmion (topological) density. Finally, there is 
a competition between Coulomb energy, which wants to 
make skyrmions of infinite size, and the Zeeman energy, 
which wants to shrink the skyrmions to a point; this de- 
termines the size of the skyrmion. (While such topologi- 
cal excitations always exist, they are not necessarily the 
lowest energy charged excitations ||.) These textures 
exist also in other two-component QHE systems, such as 
quantum Hall bilayers, in which the layer index plays the 
role of a pseudospin ||. 

Here we study skyrmions in multi-component QH sys- 
tems 0] . We consider the case where there are N degen- 
erate internal "orbitals" for each electron, leading to an 
SU(iV) symmetry. Our principal motivation for this gen- 
eralization is Si QH systems, in which, in addition to spin, 
there is a valley degeneracy for each single particle state 
H . The valley degeneracy depends on the orientation of 
the interface along which the two-dimensional electron 
gas lives. For a (100) interface, which is the orientation 
usually used in Si MOSFET QH systems, the six valleys 
are split into an upper quadruplet and a lower doublet. 
The doublet is then further split by the interfacial elec- 
tric field E, analogous to the Zeeman splitting of up and 
down spin states in a magnetic field. However, for a (110) 
interface the quadruplet lies lower and is split into two 
doublets by E, and in the case of a (111) interface all six 
valleys remain degenerate. With N v degenerate valleys 
and two spin states, the appropriate symmetry group, 
in the absence of Zeeman and other symmetry-breaking 
terms, is SU(iV) where N = 2N V . 

Our main results are as follows. We derive the en- 



ergy functional for SU(Af) QH skyrmions at filling factor 
v = 1, 2, . . . , N — 1 and show that the energy required to 
create a skyrmion-antiskyrmion pair is half the energy to 
create a (polarized) quasielectron-quasihole pair. In ad- 
dition to their charge, the skyrmions are characterized by 
N — 1 quantum numbers. For Si (100) QH systems for 
example, this implies that for small enough symmetry- 
breaking terms SU(4) skyrmions would be the lowest en- 
ergy charged excitations at v = 1,2,3. However, this 
limit is hard to reach. Better candidates for SU(4) and 
SU(6) skyrmions are Si (110) and Si (111) QH systems 
respectively. Multilayer QH systems offers another pos- 
sibility to realize N > 2 skyrmions. We consider briefly 
symmetry breaking terms and determine the energies for 
valley SU(2) skyrmions as a function of the strength of 
an easy axis symmetry breaking term. 

SU(N) QH Skyrmions - We first generalize the mi- 
croscopic derivation of the skyrmion energy and charge 
density of Yang et al. || to the case of N electron "fla- 
vors" . This will allow us to map the low-energy dynamics 
onto a nonlinear sigma model. Let fi — 1, . . . , N<p index 
the degenerate lowest Landau level single particle states 
on the plane; N,p is the total number of flux quanta. 
Consider now the Slater determinant state | [0] ) = 
rC=i 11^=1 i>ta I ) 01 integer filling factor v = m < N. 
We define the /irst-quantized operators S a p(i) in terms 
of their action in a flavor-diagonal basis for a given elec- 
tron i: S a f}(i) annihilates any state in which the flavor of 
electron i is not /? and otherwise changes the flavor from 
[3 to a. Here i € {1, .. . , A/"}, where Af = mJV^ is the 
total electron number. These operators satisfy SU(A) 
commutation relations [S a p, S p \] = S a \S p p — S p pS a \ for 
each electron and commute for different electrons. Note 
that is invariant under the U(m)x U(A — m) sub- 
group generated by S a b and S a i, where a, b 6 {1, . . . , m}, 
a, b € {m + 1, . . . , N}, and S a p — J2i S a p(i)- Following 
and H, we define the flavor texture 



*„[?] ) = exp Jd 2 r (q a - a {r) S- aa {r) 



q*a- a (r)S a - a (r)) |* m [0]) 
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where q a a(f) is a complex m x (N — m) matrix- valued 
function of r and q a a(k) is its Fourier transform, 



(2) 



and O denotes the lowest Landau level projection of O. 
The projected density operator satisfies 



exp(-ifc • ri) = Pi(k) 

p t (k) p t (p) = exp(|fe • pf + | • k hpt^p^k + p) (3) 



with I = y/hc/eB the magnetic length and k A p = z 
k x p. The ground state is homogeneous, hence 



p(fe)i* m [o]) =<j(fcyvi* m [o]> , 



(4) 



where p(k) — J2i Pi(k) is the total projected density op- 
erator. The q aa are complex coordinates parametrizing 
the coset space U(7V)/[U(m)xU(A — m)], which is the 
target space of the sigma model. Since n 2 of the target 
space is Z, topological excitations - skyrmions - charac- 
terised by one integer topological charge exist for any N 
and m. For m = 1 the target space is CP^ -1 , which for 
N = 2 is S 2 . 

Following [0] we calculate the excess charge density 
and the energy in the flavor texture state in a gradient 
expansion of q: 



5p(r) = ( # m [g] | p(r) | * m [q] ) - ( * m [0] | p(r) | * m [0] ] 



= — V 9aa A Vg* 5 + . . . 
E= (V m [q]\V\y m [q]) - (tf m [0]|F|tf m [0]> 



where 



2p° s J d 2 r \Vq a 
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v(k)p(-k)m 



dk k 3 h(k)v{k) 



(5) 
(6) 

(7) 
(8) 



h(k) = ^ (y m [o]\p(-k)m\y m [o]) . 

Here v(q) is the Fourier transform of the interparticle 
potential. For Coulomb interactions p^ — e 2 /(16\/27re^). 

To go beyond the gradient expansion is a tedious 
proposition. However, we may recognize the above re- 
sults as those obtained from expanding the energy E and 
topological density J° of the U(JV)/[U(m)xU(iV - m)] 
sigma model. A general element in the coset space may 



be written as R = U^AU, where U = exp ( ^ | ^ j and 



dH 1 "™ -i°-J" 



, i i is the matrix q aa and and 1„ is 

(1) \ U — liV-m / 

the p-dimensional unit matrix). The topological current 
is written 



J» = T7r ^ lJX Tv(Rd u Rd\R) 

l07T 



(9) 



the total topological charge, Q — Jd^J , is an integer. 
The energy functional is given by 



E Q [R] = \pl Jd 2 rTr{VRf 



(10) 



Our results are consistent with these, upon identifying 
5p = J°, to lowest order in the gradient expansion. As 
shown by ||, the long-ranged Coulomb interaction ap- 
pears at higher order, and a more accurate energy func- 
tional is therefore 

E[R] = \p° a Jd 2 rTr(VR) 2 

+ i fd\ fd 2 r' J°(r)v(r -r')J°(r') . (11) 



The dynamics are determined by the kinetic contribu- 
tion 

o ' c 

(* m [«]|i^l*mM) =^ J d 2 rq* a - a d t q a - a + ... (12) 

~^pjd 2 r J duTi(Rd u Rd t R) , 

where R(u = 0, t) = A and R(u = 1, t) = R{t). The last 
line of the above equation is the appropriate generaliza- 
tion of the m = 1 nonlinear sigma model kinetic term 
(see @). 

Skyrmion and quasiparticle energies - We now consider 
the energies of skyrmions. From Tr \d^R±i e^Rd u R\ 2 > 
we find Eq > 4717?° |Q|. The minimum energy for a 
skyrmion (or antiskyrmion) with charge Q = ±1 (electric 
charge ±e) at ^ = m is thus i? s k = 47rp° = iy / ^(e 2 /e^) 
(assuming Coulomb interactions). Comparing this to the 
energies of a polarized quasielectron, E c — 0, and a po- 
larized quasihole, E^ = y / ^(e 2 /e£), we see that, as in the 
SU(2) case, the gap to creating a skyrmion-antiskyrmion 
pair is half the gap to creating a pair of polarized quasi- 
particles. Note that the energies are independent of N 
and of m. We conclude that, for small enough symme- 
try breaking terms, skyrmions will be the lowest energy 
quasiparticles for v = m = 1, 2, 3, . . . , N — 1. Without 
the Coulomb term, the sigma model is scale invariant 
and skyrmions of any size with energy E — 47rp°|Q| ex- 
ist. Coulomb interactions make the skyrmions infinitely 
large (the energy being the same). Symmetry breaking 
terms imply a cost to texturing spins, hence they favor 
small skyrmions. The energy is then larger than the lower 



2 



bound and the size of the skyrmion is determined by a 
competition between the Coulomb interaction and the 
symmetry breaking terms. 

In addition to charge, the SU(2) skyrmions are charac- 
terized by their spin, i.e., by their S z eigenvalue. For 
SU(iV) there are N — 1 commuting generators, hence 
these skyrmions are characterized by N — 1 quantum 
numbers. Explicit symmetry breaking terms may of 
course reduce the number of conserved quantities. 

Experimental Realizations - We have shown that 
SU(iV) skyrmions are the lowest energy quasiparticles 
in any iV-component QH system provided the symme- 
try breaking terms are small enough. In the ordinary 
Si MOSFET QH systems, based on the (100) interface, 
the combination of spin and valley degeneracy leads, as 
pointed out above, to an SU(4) symmetry. However, it 
is important to note that the Zeeman term that breaks 
the SU(2) sp i n subgroup under normal conditions is too 
large for spin skyrmions to be the lowest energy charged 
excitations. If g = gfi B B/(e 2 /el) > 0.054, then the low- 
est energy quasiparticles are fully polarized 0]. In GaAs, 
this criterion is not met until B w 25T, but in Si, where 
g ~ 2, the critical field is much lower: B c sa 1.5T. Thus, 
under normal conditions, the quasiparticles are fully spin 
polarized throughout the quantum Hall regime and one 
may not be able to reach the SU(4) skyrmion regime. 

We believe the chances are best to realize N > 2 
skyrmions in QH systems based on the (110) or (111) 
interfaces of Si. It is, for example, possible to grow high 
quality SiGe heterostructures of this type, and the sym- 
metry breaking terms are believed to be small and tun- 
able. This would give four and six degenerate valleys re- 
spectively and hence lead to SU(4) and SU(6) skyrmions. 
Multilayer QH systems, possibly with spin, offer another 
possibility to realize N > 2 skyrmions, if the symmetry 
breaking terms can be made small enough. 

What happens in a particular system depends crucially 
on the type and size of the symmetry breaking terms. 
As an example of a possible scenario we discuss in some 
detail the (100) Si case, even though it may be hard to 
reach the SU(4) skyrmion regime in such systems. 

SU(4) and (100) Si QH Systems - Ignoring symmetry 
breaking terms there are four degenerate lowest Landau 
levels and skyrmions are the lowest energy charged exci- 
tations at v = 1,2, 3. If we denote the valley index by 
left (L) and right (R) and spin by up (f) and down (j), 
a basis of states is { | L \ ) , | L J. ) , | R | ) , | R J ) }, and 
the three commuting SU(4) generators can be choosen 
as <r zR = diag(0,0, 1, -1), g z l = diag(l, -1, 0, 0) and 
Tj = diag(— 1, 0, 1, 0). (y Z R/L measure the z-component 
of the spin in valley R/L, tj is the difference in the num- 
ber of spin | electrons in the R and L valleys. 

We assume the most important symmetry break- 
ing terms in Si are: the spin Zeeman term, H g = 
\g^BB(a Z B. — a z i), which affects only the SU(2) sp i n sub- 
group (and breaks it to U(l)) and a term H w , which 
affects only the SU(2) va iicy subgroup and breaks it to 
a U(l) easy axis symmetry |1C|] . (A possible choice is 



H w = w X^( r Tfc+ r l<0( T tfc+i+' r l*:+i)' w <Q and k num- 
bers the orbital states in the lowest Landau level. An- 
other possibility is (|l5|) below.) The Zeeman term causes 
the two spin { lowest Landau levels to have lower energy 
than the two spin ] ones. Since H w provides an easy axis 
anisotropy, the L/R valleys have the same energy. These 
symmetry breaking terms commute with the three diag- 
onal generators <tl, ctr, tj, which thus still give good 
quantum numbers. 

Qualitatively, the behavior of the groundstate and the 
excitations as functions of the symmetry breaking param- 
eters g and w is as follows. For v = 1, the groundstate is 
two-fold degenerate: either all L | or all R J, states are 
filled. At v = 2, all L J. and R [ states are filled and at 
v = 3, L t or R t are filled in addition. In each case, for 
small enough g and —w skyrmions are the lowest energy 
charged excitations. However, for large enough g and/or 
— w the skyrmion creation energy exceeds the energy to 
create polarized quasiparticles. Using Hartree-Fock, as 
for the SU(2) case [ fTl) , it should be possible to obtain 
the energies and quantum numbers for the skyrmions as 
functions of g, w and in particular to determine the range 
of parameters for which the skyrmions are the lowest en- 
ergy charged excitations. For large —w, the problem re- 
duces to the spin SU(2) case whereas for large g the 
problem instead reduces to the valley SU(2) case which 
we now discuss. 

Silicon Valley SU(2) - If due to the largeness of g the 
spins are fully polarized, valley-textured quasiparticles 
may be the lowest energy charged excitations in Si MOS- 
FET quantum Hall systems. For Si (100) interfaces, the 
interface potential v(z) leads to a splitting Jl2| 



AE = 



(13) 



of the lower doublet, where a ~ 0.23 A, and where the av- 
erage (dv/dz) is computed with respect to the envelope 
function A (z), as described in ref. fl^j . This is analo- 
gous to the usual Zeeman splitting, and the energy gap 
as a function of AE has been calculated previously [[[J. 
Without recourse to techniques such as nuclear magnetic 
resonance, which has been used to detect spin skyrmions, 
it may be only through the energy gap, as measured in 
transport experiments through the activated behavior of 
p XXl that valley skyrmions can be detected. Varying the 
interfacial potential v(z) with an applied electric field or 
by pressure will lead to a change in AE and hence of the 
skyrmion-antiskyrmion creation energy. 

Typically, though, the Si valley splitting AE is quite 
small - on the order of about 2 K - about 1% of the 
Coulomb energy e 2 /e£ in Si at B = 10. It seems 
likely then that disorder-induced intervalley scattering 
will dominate this Zeeman term. This is described by a 
Hamiltonian of the form 



JVi, 



>,)] 



(14) 
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where are the valley raising and lowering operators. 
The disorder thus enters as a random in-plane magnetic 
field which breaks SU(2) va ii ey down to an easy axis U(l) 
jl0|,[l3|]. As we are primarily interested in the effects of 
this internal (rather than translational) symmetry break- 
ing, we simplify the model by annealing the disorder in 
both space and (imaginary) time, and consider the follow- 
ing phenomenological Hamiltonian, which contains the 
correct symmetry breaking terms, 

H = \ [d\[dV V(r - r') : (p(r) - p ) (p(r) - p a ) : 



+\ d 2 r d 2 r'W(r-r'):T z (r)T z (r'): , (15) 



where 

E <(»")^ a (»")«<a^a- (16) 

Here we suppress the spin label (we assume full spin 
polarization); a = ±1 is the valley label, and <p^{r) is 
the normalized symmetric gauge wavefunction of angu- 
lar momentum fi in the lowest Landau level. We take 
V(r) = e 2 /er and 



W(r - r') = [W Q r 2 S(r - r') + W x V 2 8(r - 



e 2 

r')] C Tt 



(17) 



Wq and W\ are phenomenological symmetry-breaking 
parameters. In the easy axis case, W\ < 0, the symmetry 
breaking term leads to a gap — -Wi^j in the dispersion 
relation for the valley waves . Wq is a local term that 
does not affect the form of the quasiparticles. 

We now solve for the Hartree-Fock valley-textured 
quasiparticles, following the work of Fertig et al. [ [Tl|| . 
The Hartree-Fock wave functions for the skyrmion and 
anti-skyrmion are: 

oo 

|sk> = l[Ml + + v^l +h _)4_\0) 



p=0 



I sk) = Y[Ml + + v^l_ h _)\o) , 



(18) 



where \u^\ 2 + \v^\ 2 = 1 and u ll ,v il are obtained by solv- 
ing the Hartree-Fock equations self-consistently. Setting 
Ufj, = 1, we obtain the conventional quasiparticle and 
quasihole states | ^± ) , with energies 



2vr e£ 



E 4 
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7T 



e 

Ti ' 



(19) 
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FIG. 1. Skyrmion (lower) and antiskyrmion (upper) ener- 
gies versus Wi for three values of Wo- The solid and dashed 
lines correspond to the polarized quasiparticle energies of eq. 

E 

relative to the energy of the polarized ground state 
| ^[o] ) = J]^° =0 | 0) . Our Hartree-Fock results are 
plotted in figures 1 and 2. We find that valley skyrmions 
have lower energy than the polarized quasiparticles if 
-0.117 ~ Wic < Wi < 0. (This is true for any W 
as expected since this is a local term ||.) 

Discussion - To summarize, TV component QH systems 
have an SU(7V) symmetry and SXJ(N) skyrmions are the 
lowest energy charged excitations at v = 1, 2, . . . , N— 1 in 
the symmetric limit. N > 2 skyrmions may exist as val- 
ley skyrmions in (110) and (111) Si systems, as combined 
spin and valley skyrmions in standard (100) Si systems 
or in multi-layer QH systems. The most promising can- 
didate probably being (110) or (111) SiGe systems. 

What happens in a particular system depends crucially 
on the type and size of the symmetry breaking terms and 
experimental signatures will depend on the coupling to 
external probes. The energy gap to creating the lowest 
energy charged quasiparticles can be measured in acti- 
vated transport, and for SU(2) sp j n skyrmions the spin 
can be obtained from such measurements by changing 
the in plane magnetic field. An external interfacial elec- 
tric field plays a similar role for SU(2) va n C y skyrmions 
in Si, although in this case the energy depends on the 
field in a more complicated way. This can be used as 
a signal for SU(2) va n cy skyrmions and, in combination 
with an in plane magnetic field, to study crossovers be- 
tween skyrmions with various quantum numbers in a 
limit where g, —W\ and the valley Zeeman term are all 
small. 
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